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One of the simplest quantum spin systems is the 
spin-i XY chain introduced by Lieb et al [1], who 
pointed out the relation between the spin model and 
noninteracting spinless fermions that makes many 
properties of the spin model amenable to analytic cal- 
culations. Dynamic spin correlations are of particular 
interest since they can be measured by neutron scat- 
tering or magnetic resonance techniques in a number 
of materials which are reasonably well described by 
different special cases of the spin-i XY chain (e.g., 
CSH2PO4, PbHP04, PrClg, CsCuClg etc). Unfortu- 
nately, the calculation of dynamic spin correlations 
can be nontrivial even in the simple case of the XY 
chain, and the complete wavevector and frequency de- 
pendences of the dynamic structure factors or the dy- 
namic susceptibilities have been obtained only recently 
using a numerical approach [2]. In what follows we 
present the results for the time-dependent two-spin 
correlation functions of the spin-i isotropic XY chain 
in a transverse field obtained numerically and apply 
these findings to discuss the temperature dependence 
of the spin-spin relaxation time in PrCls [3] which was 
not known heretofore due to lack of information about 
the correlation functions for the spin-^ isotropic XY 
chain. 

We consider N spins one-half governed by the fol- 
lowing Hamiltonian 
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if = f7 V + J V («+i + 44+1) 



where n is the transverse field and J is the ex- 
change interaction. Using the numerical approach ex- 
plained in detail in Ref. [2] we compute the time- 
dependent two-spin correlation functions (•s?(i)s^i „), 
a/3 = XX, xy, zz. Here the angular brackets stand for 
the canonical thermodynamic average with the inverse 
temperature (3 = l/ksT . In figure A we show the nu- 
merical results for the time dependence of the real part 
of the XX autocorrelation function for fi = 0, J = — 1 
at different temperatures. In our numerical compu- 
tations we take N = 400, j = 51. Downward and 
upward triangles, circles, and squares correspond to 
/3 = 5, 1, 0.1, and 0.00001, respectively We also 
plot the exact analytical result for infinite tempera- 
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ture /3 = [4] by a solid curve. Dot-dashed, short- 
dashed, long-dashed, and dotted curves correspond 
to the asymptotic result [5] for (3 = 5, 1, 0.1, and 
0.00001, respectively. From figure A one can see an 
excellent agreement between the numerical and ana- 
lytical results (only the slopes of asymptotics should be 
compared with numerical results). At finite tempera- 
tures (sj(i)sj) decays exponentially (the logarithms of 
the asymptotic results have the same slopes as the log- 
arithms of the computed correlation functions); with 
increasing temperature one observes Gaussian decay 
over a time interval of increasing length (the plotted 
data for (3 = 0.00001 show that the correlation func- 
tion at this temperature does not reach its asymptotic 
regime within the time range displayed in the figure 
A). ^ 

The time-dependent two-spin correlation functions 
yield the dynamic structure factors or the dynamic 
susceptibilities that can be measured in materials 
which are reasonably well described by the spin-i XY 
chain. As an application of our study on the time de- 
pendence of the spin correlations we reconsider the 
theoretical prediction for the temperature dependence 
of the spin-spin relaxation time in PrCls [3] . The Pr- 
Pr interaction in this compound as derived from mea- 
surements of static quantities can be reasonably well 
described by the model considered here with the ex- 
change interaction J/Ub — 2.85 K and fi = 0. The 
spin-spin relaxation time T2 is related to the xx time- 
dependent autocorrelation function 



1 
T'2 



dt(s?(i)4> 



In figure B we show the result of our long-chain calcu- 
lation of the dependence of T2 on temperature. T2 is 
an increasing function of temperature and crosses over 
to a constant at about 1 K. The plateau in the tem- 
perature dependence of T2 can be understood from the 




kT/J 



results presented in figure A. For the time range dur- 
ing which (s|(t)s|) contributes appreciably to the inte- 
gral, it does not depend on temperature for sufficiently 
high temperatures. As one can see in figure A the data 
for (3 < 1 coincide in the region where {sj{t)s^) has 
non-negligible values; the diS'erences between high and 
infinite temperatures occur only at long times, where 
{sj{t)s^) is already very small. 

Qualitatively, our result for the temperature depen- 
dence is similar to that of an earlier short-chain calcu- 
lation [3] and to the behaviour observed experimen- 
tally [3]. However, there are significant discrepan- 
cies in the value of the crossover temperature. The 
short-chain calculations yields a crossover at less than 
2.5 K; additional approximations (coarse-graining in 
the frequency uj) were necessary to extrapolate from 
A^ = 10, 12 to A^ = cxD [3]. The experimental crossover 
temperature of 6 K is still much higher. This indicates 
(as already suspected earlier [3]) that while the static 
thermal behaviour of PrCla may be describable by the 
isotropic XY chain with J/fc^ — 2.85 K, the dynamics 
is probably more complicated. 
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